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1. Analysis of the Continuous-Data Instrument Pointing System 

The objective of this chapter is to investigate the performance of 
the simplified continuous-data model of the Instrument Pointing System 
(IPS). Although the ultimate objective is to study the digital model 
of the system, knowledge on the performance of the continuous-data model 
is important in the sense that the characteristics of the digital system 
should approach those of the continuous-data system as the sampling period 
approaches zero. 

The planar equations which describe the motion of the Spacelab using 

inside-Out Gimbal (iOG)- for the pointing base are differentia] equations 

of the fourteenth order. A total of seven degrees of freedom are 

represented by these equations. Three of these degrees of freedom (X , Z , 

s ’ s ’ 

$ ) belong to the orbitor, and three degrees of freedom (X , Z , 4') are 

m m 

for the mount, and the scientific instrument (Sl) has one degree of free- 
dom in e. 

A simplified mode] of the IPS model is obtained by assuming that all 
but motion about two of the seven degrees of freedom axes are negligible. 
The. simplified IPS model consists of only the motion about the scientific 
instrument axis, £, and the mount rotation 'F. The block diagram of the 
simplified linear IPS control system is shown in Figure 1-1. 

in this chapter we shall investigate only the performance of the 
continuous-data iPS control system. The signal flow graph of the 
continuous-data IPS control system in Figure 1-1 is shown in Figure 1-2. 

The characteristic equation of the continuous-data IPS is obtained 
from Figure 1-2. 


A = (1 - K^K^)s^ + (K^ + - K,K^K^)s^ + (K^ - + K^K^K^)s^ 
























0 


+ (K,K^ - K|K^K^ + K|KjKj + V^K^) 5^ + (K, + K^K^K^js ^ = 


( 1 - 1 ) 

The system parameters are: 

Kq = 8 X 10 ^ n-m 

. 4 

= 6 X ]0 n-m/sec 

= 0.0012528 
= 0.0036846 
Kj = var i abl e 
= 0.80059 
= 10798.49 
= 1.1661 
= 0.0000926 

Equation ( 1 - 1 ) is simplified to 

s^ + 16.704s^ + 222. 63s^ + (1.706 -t- 27. 816 10’'^K|)s^ 

+ (4.11 + 174.7 10'\|)s + 513.855 10‘^Kj = 0 ■ (1-2) 

The root locus plot of Eq. ( 1 - 2 ) .when K| varies is shown in Figure I-3. 
It is of interest to notice that two of the root loci of the fifth-order 
IPS control system are very close to the origin in the left-half of the 
s-plane, and these two loci are very insensitive to the variation of Kj . 

The characteristic equation roots for various values of Kj are tabulated 
in the fol lowing: 


K. 


10 


10 ' 


10 '- 


ROOTS 

0, -0.003l4jj0. 13587, -8.3488+jl2.36l4 
-0.0125 -0.00313984HJ0. 135873, -8.34260+^)12.3571 
-0.1262 -0.0031342+j0. 135879, -8.28577+j]2.3.192 
“1.38146 -0.0031398+jo. 135882, -7.658I 3+j 1 1 .9457 




ane I 






















K 


ROOTS 


2X10^ 

-3.08107, 

-0.00313629+30.135883, 

-6.80833+311.5845 


5x10^ 

-9.07092, 

- 0.0031401 IjJO. 135881 , 

- 3 . 81340 +jll .7806 


2 x 10 ^ 

-19.7201 , 

-0. 003 140254J0. 135880, 

-1 .511 184 jl 6.7280 


107 

- 14 . 5468 , 

- 0 . 00314020 + 30 . 135881 , 

- 1 .07547+313.7862 


5x10^ 

-27.2547, 

- 0 . 00314063 + 30 . 135889, 

5 . 27850 +J 21 .963 


10^ 

-34.-97, 

-0.00314035+30. 135882 , 

8 . 69964 jJ 27.2045 


The continuous-data IPS control system is 

asymptotically stable 

for K| 

1 ess than 1 . 5 x 10 ^. 



It 

is of interest 

to investigate the response of the IPS control 

system 

due to its own 

initial condition. The 

transfer function 

between 

and 

its initial condition e for K, = 5 x 10 ° 

0 1 

i s 


e(s) _ 

s'(s 3 

+ 16 . 7 s^ - 893 s - 5581 . 

77 ) 

( 1 - 3 ) 


+ 16 . 7 s^ 

+ 222 . 63 s^ + 1392.55^ + 

12 . 845 s + 25.693 


where we have considered that e(O) ~ is a unit~step function input 
applied at t = 0; i,e., e^(s) = 1/s. 

it is interesting to note that the response of B due to its own 
initial condition is overwhelmingly governed by the poles near the 
origin. In this case, the transfer function has zeros at s = 8O.5, 
“I2.375+J23.75. The zero at s= 80.5 causes the response of e to go 
negative first before eventually reaching zero in the steady^-state, 
as shown in Figure 1 -^. The first overshoot is due to the complex 
poles at s = -3 . 8 l 3 ^jJ 1 1 .7806. Therefore, the ei genvalues of the 
closed-loop system at s = -3 . 8 l 3 ^jJ 1 1 *7806 controls only the transient 
response of s(t) near t = 0, and the time response of e(t) has a long 
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period of oscillation. The eigenvalues at s = -0.00313 + j0.1358 due to 

the isolator dynamics give rise to an oscillation which takes several 

minutes to damp out. The conditional frequency is 0.1358 rad/sec, so the 

period of oscillation is approximately seconds. Figure 1-5 shows the 

response of e(t) for = 1 rad on a different time scale over 100 seconds. 

Although the initial condition of 1 radian far exceeds the limitation under 

which the linear approximation of the system model is valid, however, for 

linear analysis, the response will have exactly the same characteristics 

but with proportionally smaller amplitude if the value of e is reduced. 

o 

Since the transient response of the system is dominated by the 
isolator dynamics with eigenvalues very close to the origin, changing 
the value of Kj within the stability bounds would only affect the time 
response for the first second as shown in Figure 1-4, the oscillatory 
and slow decay characteristics of the response would not be affected in 
any significant way. 
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2. The Simplified Linear Digital tPS Control System 

In this chapter the model of the simplified digital IPS control 
system will be described, and the dynamic performance of the system will 
be analyzed. 

The block diagram of the digital IPS is shown in Figure 2-1, where 
the element S/H represents sample-and-hold . The system parameters are 
identical to those defined in Chapter 1. An^ equivalent signal flow graph 
of the block diagram in Figure 2-1 is drawn as shown in Figure 2-2. 

Applying Mason's gain formula to Figure 2-2 with 0 (s) and fi.(s) as 

A A 

outputs, 0 ^*(s) and fl^*(s) as inputs, we have 


G,(s) 


Gh(s) 

As^ 


KiT 


KJ 


where 


A^ = 1 + K^s + K^s 


-2 


S(s) “ 


1 - e 


^ 2 - 

-Ts 


h 

^*^7 As 

" Ka) 0 aMs) 

( 2 - 1 ) 

G (s) 

•)K (A 

‘ As 

- Kij)0A*(") 

( 2 - 2 ) 


(2-3) 


-2 

3 S 

(2-4) 



(2-5) 



Lett i ng 




“ e [(1 - K^)s^ + K^s + K^) 

3] 


= ,(s) = - K,) = ~ (2-6, 


^ S^((l .- K^Kg)s^ + K^s + kJ 


G 2 (s) 


Gj.^(s) G.(s) 


(2-7) 


and taking the z-transform on both sides of Eqs.(2-l) and (2-2), we have 


KjT 


J^a(2) KjGj (z)<;^a(z) - (f^o 


( 2 - 8 ) 


K.T 


0a(2) = - (Kq +^)G2(z)0 (z) 


(2-9) 
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From these two equations the characteristic equation of the digital system 


is found to be 


A(z) = 1 + (z) + (Kg + ° 


( 2 - 10 ) 


where 


Gl(z) = ^(G^(s)] = K^(l - 


(1 - K^)s + K^s + 

s^[(l - K^Kg)s^ + K^s + K ] 


( 2 - 11 ) 


G^iz) = ^( g ^( s )/ s ) 


( 2 - 12 ) 


The characteristic equation of the digital IPS control system is of 
the fifth order. However, the values of the system parameters are such 
that two of the characteristic equation roots are very close to the z = 1 
point in the z-plane. These two roots ^are rtis i de the unit circle and they 
are relatively insensitive to the va.lues of K| and T, so long as T is not 
very large. The sampling period appears in terms such as e"®-009^T 
is approximately one unless T is very large. 

Since the values of K 2 and are relatively small, 


G^ (z) - K^d - 


IJS 1 ViljV T 


(1 - K^Kg)s^ 


- SH ^ ’ 


= 2.79x10 


z - 1 


( 2 - 13 ) 




K K (1 - z = 2 . 79 x 10 “^ ^ (2-U) 

’ 2(z - 1)2 


Substituting G^(z) and G 2 (z) into Eq. (2-10), the characteristic equation 
of the digital IPS is approximated by the following third-order equation: 


, fK.K T 


K K T^ K.K T' 


K,K T - 3lz^ + [■ 
1 P [ 


K K.T- 
P I 


- 2K.K T + 3 z 
1 P 


(2.-15) 



where = 2.79^10 , and it is understood that two other characteristic 
equation roots are at 2 = 1 , It can be shown that in the limit as T 
approaches zero, the three roots of Eq . (2-15) approaches to the roots of 
the characteristic equation of the continuous-data IPS control system, and 
the two roots near 2 = 1 also approach to near 5 = 0, as shown by the 
root locus diagram in Figure I- 3 . 

Substituting the values of the system parameters into Eq. (2-15), we 

have 

A(z) = + (III.ST^ + 16.74T - 3)z^ + (3 - 33. ^8 t + 1 . 335><1 T^) z 

+ (1 .395xl0“^K|T^ + 16.74T - lll.6T^ - l) = 0 .(2-16) 


Applying Jury’s test on stability to the last equation', we have 


or 


K K.T^ > 0 
P I 


K| > 0, since T > 0. 


(1) A(l) > 0 
Thus, 

(2) A(-l) <0 or -8 + 66 . 96 T < 0 

Thus, T < 0.12 sec, 

(3) Also, |aQ|<a^, or 

1 1 . 335 x 1 o'^KjT^ + 16.74T - 1!1.6T^ - ij < 1 


(2-17) 

( 2 - 18 ) 

( 2 - 19 ) 

(2-20) 


The relation between T and Kj for the satisfaction of Eq. (2-20) is 
plotted as shown in Fig. 2-3. 

(A) The last criterion that must be met for stability is 


I^qI > (2-21) 

where , 2 2 

*^0 = ^0 “ ^3 

^2 "" ^ 0®2 " ^ 1^3 

3q = 1 . 395 x 1 0"^KjT^ + 16.7^T - 111.6T^ - 1 
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= 3 - 33.^81+ 1 .395xl0“VjT^ 

32= 111.6T^ + 16.7^T - 3 
= 1 

The relation between T and K| for the satisfaction of Eq. (2-21) is 
plotted as shown in Figure 2-3. It turns out that the inequality condition 
of Eq. (2-21) is the more stringent one for stability. Notice also that ’ 
as the sampling period T approaches zero, the maximum value of Kj for 
stability is slightly over 10^, as indicated by the root locus plot of 
the continuous-data IPS. 

For quick reference, the maximum values of for stability corres- 
ponding to various T are tabulated as follows: 

T Max. K| 

0.01 10 ^ 

0.08 7.5x10^ 

0,1 6.9x10^ 

When T = 0.05' sec, the characteristic equation is factored as 

(z - l)(z^ - 0.884z + 0.442 - 0.1744 x10”^K|) = 0 (2-22) 

Thus, there is always a root at a = 1 for all values of K| , and 
the system is not asymptotically stable. 

The root locus plot for T = 0.01, 0.08, and 0.1 second when Kj 
varies are shown in Figures 2-4, 2-5, and 2-6, respectively. The two 
roots which are near z = 1 and are not sensitive to the values of T and 
K| are also included In these plots. The root locations on the root loci 
are tabulated as follows: 
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T = 0.01 sec 


•^1 

• 

ROOTS 

0 

1.00000 

0.91072 + jo.l 19788 

10^ 

0.986'i 

0.917509 + jo. 115822 

5x10^ 

0.907357 

0.957041 + jO.l 14946 

lo7 

0.853393 

0.984024 + jO. 137023 

5x10^ 

0,73^201 

1 .04362 + jo. 224901 ■ 

10^ 

0.672306 

1 .07457 + jo. 282100 

T = 0.08 



ROOTS 

0 

0.999972 

-0.0267061 + jO. 611798 

10^ 

0.888688 

0.0289362 + j 0,583743 

5x10^ 

0.0273573 

0.459601 +jO. 665086 

7.5x10^ • 

-0.156092 

0.551326 + jo. 851724 

10^ 

-0.253228 

0,599894 + jO. 989842 

108 

-0.770222 

0.858391 + j2. 83716 



E:UG£NE DJETZCEN CD 











































ROOTS 


0 

0.990937 

10^ 

0 . 860648 

5>^10^ 

-0.417697 

7x10^ 

-0.526164 

107 

-0,613794 

10« 

-0.922282 


“0.390^69 + jo. 522871 
"0. 325324 + jO. 495624 
0.313849 + jO. 716370 
0.368082 + j 0.938274 
0.411897 + ji . 17600 
0.566141 + j3. 78494 


The linear digital IPS control system shown by the block diagram of 
Figure 2-1, with the system parameters specified in Chapter 1, was 
simulated on a digital computer. The time response of e(t) when the- 
initial value of e(t), is one, is obtained. Figure 2-7 illustrates 
the responses of e(t) for T = 0.01, Kj = 5=<10^, and T = 0.1, Kj = 10^ 
and K| = 10^. Similar to the continuous-data IPS control system analyzed 
in Chapter 1, the time response of the digital IPS is controlled by the 
closed-loop poles which are very near the 2 = 1 point in the z-plane. 
Therefore, when the sampling period T and the value of Kj change within 
the stable limits, the system response will again be characterized by 
the long time in reaching zero as time approaches infinity. 

The results show that as far at the linear simplified model is 
concerned, the sampling period can be as large as 0.1 second, and the 
digital IPS system is still stable for less than 10^. 
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3. Analysis of Continuous-Data IPS Control System With Wire Cable 
Torque Nonlinearity 

in this chapter we will conduct a stability analysis of the simplified 
IPS control system with the nonlinear characteristics of the torques caused 
by the combined effect of the flex pivot of the gimbal' and by wire cables, 
for transmitting power, etc., across the gimbal to the experiment, as shown - 
i n Figure 3~1 • 

The flex pivot torque disturbance has been modeled^ as a linear spring 
with slope Kpp, as shown in Figure 3“2. The wire torque disturbance is' 
modeled as a nonlinear spring, with a slope of Ky.j., as shown in Figure 3-3. 
The total torque disturbance associated with the flex pivot and wire cables 
is summed up as shown in Figure 3"^- The continuous-data IPS control system 
with the nonlinear torque characteristics Is shown in Figure 3~5. 

The objective of the analysis in this chapter is to study the condition 
of sustained oscillation of the nonlinear IPS control system shown in 
Figure 3-5. 

The Describing Function of the Wire-Cable Nonlinearity 

it was pointed out in the above discussion that the wire cable 
nonlinearity can be modeled by either the arrangement shown in Figure 3-^(a) 
or Figure 3-^(b). For the model of Figure 3“^{3), a relay characteristic 
is present between e and T , and, in addition, a linear gain of K, 
exists between e and T , where K denotes the gain constant due to the flex 
pivot torque. Using this model, one can conduct a describing function 
analysis with the relay nonlinearity, but the linear system model is altered 
by the addition of the branch with the gain of + Kpp. However, careful 
examination of the block diagram of Figure 3~5 reveals that the branch with 






Figure 3-3. Wi re 


/ 


Figure 



(b) 


k. Combined flex pivot and wir 
characteri sties. 























the gain of is parallel to the branch with the gain which has 

a magnitude of 8xl0^ N-m/rad. Since the value of K,,^ lies between 0.25 and 
25 N-m/rad, and the maximum value of K^p is In the order of several hundred, 
it is apparent that the value of will be predominant on the system 
performance. This means that the linear transfer function will, not change 
appreciably by the variation of the .values of and Kpp. 

Prediction of Self-Sustained Oscillations With the Describing 
Function Method 


From Figure 3"5 the equivalent characteristic equation of the nonlinear 
system is written as 

1 + N(£)G^^(s) = 0 (3-1) 

where N(e) is the describing function of the relay characteristic shown in 
Figure 3~^(a), and is given by 


N(e) = ^H^.p/Tre 

The transfer function G^^(s) is derived from Figure 3~5, 
R = 0.00139^6(5^ + 0.0012 528s^ + 0.00368^6s^) 

eq^^^ A(s) 


where 

A(s) = + I6.704s^ + 222. 6s^ + (2. 78 ] xl Kj+l . 706 ) 

+ (^{.11 + 1 .7^7xl0"^K|)s + 5]3.855 x1o'^Kj 
A necessary condition of self-sustained oscillation is 


(3-2) 


( 3 - 3 ) 


-(3-4). 


G (s) = -1/N(£) 


( 3 - 5 ) 


Figure 3-6 shows the frequency response plots of G (s) for K =10^ 
67 eq I ’ 

10 , and 10 , as a function of w , and the trajectory of -1/N(e) =- 7re/4H • 

WT’ 
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the latter lies on the -180° axis for all combinations of magnitudes of e and . 

Figure 3-6 shows that for each value of Kj there are two equilibrium 
points, one stable and the other unstable. For instance, for Kj = 10^, the 
G^^(s) curve intersects the-l80° axis at w = 0.16 rad/sec and w = 0.05 rad/sec. 
The equilibrium point that corresponds to o) = 0.16 rad/sec is a stable equili- 
brium point, whereas w = 0.05 rad/sec represents an unstable equilibrium point. 
The stable solutions of the sustained oscillations for Kj = 10'’, 10^, and 10^ 
are tabulated below: 


"l 

e/H^^ (rad/sec) 

(rad) 

(arc-sec) 

m(rad/sec) 

f (sec) 
c 

10? 

7.2x10"® 

2.39x10'^ 

0.05 

0.3 

21 

10^ 

5.07x10"^ 

3. 17x10"® 

0.65 

0. 16- 

39 

lo5 

1 .13x10“^ 

8.19x10“® 

16.9 

0.14 

45 


The conclusion is that self-sustained oscillations may exist in the 
nonlinear continuous-data control system. 

Since the value of i:s very large, the effect of using various values 
of and within their normal ranges would not be noticeable. Changing 
the yalue of has a one-to-one effect on the amplitudes of oscillations 

• I 

of e and e. 

^gital Computer Simulation of the Continuous-Data Nonlinear IPS Control System 

Since the transient time duration of the IPS control system is exceedingly 
long, it is extremely time consuming and expensive to verify the self-sustained 
oscillation by digital computer simulation. Several digital computer simulation 
runs indicated that the transient response of the IPS system does not die out 
after many minutes of real time simulation. It should be pointed out that the 
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describing function solution simply gives a sufficient condition for self- 
sustained oscillations to occur. The solutions imply that there is a certain 
set of initial conditions which will induce the indicated sel f-susta ined 
oscillations. However, in general, it may be impractical to look for this 
set of initial conditions, especially if the set is very small. it is entirely 
possible that a large number of simulation runs will result in a totally 
stable situations. 

Figure 3~7 illustrates a section of the time response of e(t) from t = 612 
sec to 692 sec. The initial conditions are e(0) = lo"^ and £(0) = 10~^, and 
all other initial conditions are set to zero; K = 10^ H =1 and K. + K 
= 100. It was mentioned earlier that the system is not sensitive to the values 

\lT *^FP‘ 3-7 It is seen that the period of the oscillation 

is 1|8 seconds, or 0.13 rad/sec, which Is very close to the predicted value. 



Figure 3“7- A section of the response of s(t) 
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4- Analysis of the Digital IPS Control System With Wire Cable Torque ■ 

Nonl ineari ty 

In this chapter we will conduct a stability analysis of the simplified 
digital IPS control system with nonlinear characteristics of the torques 
caused by the combined effect of the flex pivot of the gimbal and wire cables. 

The analysis used here is the discrete describing function which will 
give sufficient conditions on self-sustained oscillations in non 1 inear d i g i tal 
systems. The block diagram of the nonlinear digital IPS system is shown in 
Figure ^f-l . For mathematical convenience, a sample-and-hold is inserted at 
the input of the nonlinear element. 

A signal flow graph of the system in Figure 4-1 is drawn in Figure 4-2. 
The z-transforms of the variables in Figure 4-2 are written. 


0^(z) = -G^(z)T^(z) 
(z)T (z) 


where 


e;(z) = -G^(z)T^(z) 


K|T 


Tc(z) - K,a^(z) + (K^ +7TT) ©A^^) 

- i\ij + Kpp)G3(z)T^(z) + N(z)e(z) 
Gjs) 


(z) 




*^7 As ^^1 ■ W 


K. 


G,(s) 


'3 


(z) = ^ 


7 As 


r G (s) 


2~(A^ - K^) 


a^U) = j[k 


G,(s) 


7 As 1 
-Ts 


- e 


= s 
A^ = 1 + K^s”^ + K^s“^ 


(4-1) 

(4-2) 

(4-3) 


(4-4) 






















36 


A = (1 ~ K^Kg) +-K2 s"' + K^s 


-2 


N(z) - discrete describing function of ideal relay 0, ~H , ) 

WT wT 

Since and are very small, approximations lead to 

. . ? . -79 [o~\ 

G,(z) . . 2 -79 10-V(z 1) 


2(z - 1)^ 

G (z) 0.000697 T^(z + 1) 




2(z - ])■ 

0.0001394T 
2 - 1 


Equations (4-1) through (h-^) lead to the sampled flow, graph of Figure 
^“3, from which we have the characteristic equation, 


A(z) = 1 + (z) + 

Equating A(z) to zero. 


K T 

K + -V 

o z-1 


^2^^^ ^\lj + G2j(z)N(z) (if-5<) 


the equivalent transfer fucntion G (z) is obtained, 

eq ’ 

— 


G^ (z) = 

1 + K,G,(z) + 


KiT 




(^-6) 


The intersect between G (z) and -1/N(z) gives the condition of self- 

eq 

sustained oscillations in the digital IPS control system. 

Figure k-h shows the G (z) plots for various values of N >-2 for K = 10^. 

eq I 

in this case, it has been assumed that the periods of oscillations are integral 

multiples of the sampling period T. Therefore, if T denotes the period of 

• c 

oscillation, T^ = NT, where N is a positive integer > 2, 

From Figure h-h it is seen that when N is very large and T is very small, 
Ggq(z) approaches G^^(s). However, for relatively large sampling periods. 
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Figure 4-3- Sampled signal flow graph of the digital IPS 
control system. 

does not approach to G (s) even for very large N. This indicates 

, 6C} 

the fact that digital simulation of the continuous-data IPS can only be 

carried out accurately by using extremel y~sma 1 1 sampling periods. 

The critical regions for -1/N(z) of an ideal relay are a family of 

cylinders in the gain-phase coordinates, as shown in Figure 4-5. These 

regions extend from -«> db to db, since the dead zone of the ideal relay 

is zero. For N = 2, the critical region is a straight line which lies on 

the -180° axis. The widest region is for N = 4, which extends from -225° 

to -135°. Therefore, any portions of the G (z) loci which do not lie in 

eq 

the critical regions will correspond to stable operations. 

It is interesting to note from Figure 4-4 that the digital IPS system 
has the tendency to oscillate at very low and very high sampling periods, 
but there is a range of sampling periods for which self-sustained 
oscillations can be completely avoided. 
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Referring to Figure k-k it is noticed that when the sampling period T 

is very small, (T ^ 10 sec approximately), the loci of G (z) for N = 2 

eq 

3, and 4 lie in their respective critical regions, and self-sustained 
oscillations characterized by these modes are possible. However, since the 
actual period of these oscillations are so small, being 2, 3, or 4 times 

I 

the sampling period which is itself less than 0.01 sec, the steady-state 

oscillations are practically impossible to observe on a digital computer 

simulation unless the print-out interval is made very small. This may 

explain why it was difficult to pick up a sel f -susta i ned oscillation in the 

digital computer simulation of the continuous-data IPS system, since a 

digital computer simulation is essentially a sampled-data analysis. 

When the sampling period is large (T ^ 9 sec approximately), the digital 

IPS system may again exhibit self-sustained oscillations, as shown by the ^ 

G (z) loci converging toward the -180° axis as T increases. However, the 

Ggq(2) loci of Figure 4-4 show that there is a midrange of T for which the 

digital system is stable. The G (z) locus for N = 2000 actually represents 

eq 

the locus for all large N. Therefore, for T = 0.1, the G (z) locus point 

eq ^ 

will be outside of the critical regions, since as N increases, the widths 
of the critical regions decrease according to 

N = even width of critical region = 2ir/N 

N = odd width of critical region + tt/N 

•Therefore, from the standpoint of avoiding self-sustained oscillations 
in the simplified digital IPS control system model, the sampling period may 
be chosen to lie approximately in the range of 0.01 to 1 second. 
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Digital Computer Simulation of the Digital Nonlinear IPS 

« 

Control System 

The digital IPS control system with wire-cable nonlinearity, as shown 
by the block diagram of Figure ^-1, has similar characteristics as the 
continuous-data system, especially when the sampling period is small. The 
digital IPS control system of Fig. ^-1, without the sample-and-hold in front 
of the nonlinear element, was simulated on the IBM 360/75 digital computer. 

With initial conditions set for c and e, typical responses showed that the 
nonlinear digital IPS system exhibited a long oscillatory transient period. 
Figure ^-6 shows the beginning portion of the response of e(t) for e(0) = lo"^, 
e(0) = lo"^, Kj = 10^, T = 0.1 sec, = 1, \<yj + K^p = 100. Figure h-J 

shows the same response over the period of 765 sec to 1015 seconds. Figure 
^-8 gives the response of e{t) for the fime duration of 1530 sec to 1725 sec,' 
and it shows that the response eventually settles to nonosci 1 latory and 
finally should be stable. 
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Figure ^-8. Response of e(t) of digital IPS control system 
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5. Gross Quantization Error Study of the Digital IPS Control System 

This chapter is devoted to the study of the effect of gross quantization 
in the linear digital IPS control system. The nonlinear characteristics of 
the torques caused by the combined effect of the flex pivot of the gimbal 
and wire cables are neglected. 

Since the quantization error has a maximum bound of +^h/2, with h as the 
quantization level, the "worst" error due to quantization in a digital system 
can be studied by replacing the quantizers in the system by an external noise 
source with a signal magnitude of jji/2. 

Figure 5“1 shows the simplified digital IPS control system with the 
.quantizers shown to be associated with the sampl e-and~hold operations. The 
quantizers in the displacement 0^, rate and torque T^, channels are denoted 
Qqj j and respectively. The quantization levels are represented by 
^o? h^, and h.^, respectively. 

Figure 5-2 shows the signal flow graph of the digital IPS system with 
the quantizers represented as operators on digital signals. Treating the 
quantizers as noise sources with constant amplitudes of +h /2, +h,/2, and 

— o ’ — I ’ 

+^h.j./2, we can predict the maximum errors In the system due to the effect 
of quantization. The following equations are written from Figure 5-2. Since 
the noise signals are constants, the z-transform relations include the factor 
'z/ (z - 1 ) . 


T (z) 
c 






= K|!i^(z) + K| 




-K^G^(s)A^ (s) 


sA(s) 


- 2 z - 

sA(s) 


K + — 
o z - 1 




+ 






s^A(s) 


s^A(s) 


'T z 


^c(-^±Tz-i 


(5-1) 

(5-2)’ 


(5-3) 




4 ?- 

a^ 






















The last two equations are written as 


where 


= g,U)[t^(z) ±^Ht) 


0^(z) = 62(2) 


T,(z) 


A, (s) = 1 + K^s + K^s 


-2 


A(s) = 1 - + K^s ’ + K^s ^ 


Gj(z) 

G,(z) 


2.78x10 \ 
z - ] 

2 .78 x10~V(z + 1 ) 

2(2 - 1 )^ 


(5-A) 

(5-5) 

(5-6) 

(5-7) 

(5-8) 

(5-9) 


Figure 5-3 gives the digital signal flow graph representation of Eqs. (5-1) 

( 5 - 4 ), and ( 5 - 5 )* Using Figure 5-3 we can analyze the worst-case errors due to 

quantization in the steady state at any point of the IPS system. 

As derived in previous chapters, the transfer functions (z) and 0 ^( 2 ) 

■< 

are given in Eqs. (5-8) and (5-9). The characteristic equation of the system 
is obtained from Figure 5~3. 


A(z) = 1 + KjGj (z) + 


K T 

K + — 
o z - 


G2(z) 


( 5 - 10 ) 


We shall now evaluate the maximum steady-state quantization errors for 

terms of the quantization levels h , h, , and h . 

' 0 I 7 

From Figure 5 - 3 , T^(z) is written 


T^(z) 


= i{+ ^ 

a(- 2 


K T 

K + — ~ 
o z - 1 




K + — 
o z - 1 


G 2 (.) 


z - 1 


(5-11) 


. The steady-state value of T^(kT) is given by the final-value theorem, 


1 im T (kT) = 1 im (1 - z"’)T (z) 




z-^l 


( 5 - 12 ) 


Substitution of Eq. ( 5 -II) into Eq . (5-12), we have 



^9 


Sitni larly, 


1 im T^(kT) = + ~ 


0a(z) 


if h 
A 2 




+ 



Then , 


1 _ h 

lim O.(kT) = lim (] - z“')0,(z) == + ^ 
k-^ ^ z-^1 A 2 


f!A<2) 



+ 


K,T 


z - 1 


+ 





(5-13) 


(5-1^) 

(5-15) 

(5-16) 


lim (kT) = lim (1 - z ^)fi (z) = 0 (5-17) 

k-4co ^ ft 


Therefore, we conclude that the maximum error in T^ due to quantization 

Is h^/2 and is not affected by the other two quantizers. The quantizer in 

th® ®A channel affects only 0^(z) In a one-to-one relation. The quantizer 

in the path does not affect either 0. , , or T . 

^ A A c 





Figure 5-3- Digital signal flow graph of the simplified digital IPS 
with quantization. 
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6. Describing Function Analysis of the Quantization Effects, of the Digital 
iPS Control System 

in this section the effects of quantization in the digital IPS control 
system are investigated with respect to self-sustained oscillations. 

Since the quantizers represent nonlinear characteristics, it is natural 
to expect that the level of quantization together with the selection of the 
sampling period may cause the system to enter into undesirable self-sustained 
osci 1 lations. 

The digital IPS control system with quantizers located in the 0 and 

A A 

channels is shown in Figure 5~1 • We shall consider the effects of only one 
quantizer at a time, since the discrete describing function method is used. 

With reference to the signal flow graph of Figure 5~2, which contains ail 
the quant i zer.s ,• we can find the equivalent characteristic equation of the 

1 

system when each quantizer is operating alone. Then, the equivalent linear 
transfer function that each quantizer sees is derived for use in the discrete 
describing function ana.lysis. 

Quantizer in the 0„ Channel 

Let the quantizer fn the 0^ channel be denoted by Q^, as shown in Figure 
5“2, and neglect the effects of the other quantizers. Let the discrete 
describing function of Q^ be denoted by Qq(z). From Figure 5-2 the following 
equations are written: 


TcU) = 


K T ^ 

K + ^ ' 

o z - 


Q^(z)0^(z) + K,fi^(z) 




rK^Gp^(s)A^(s). ^ 


s^A 


s^A 


T^U) 


( 6 - 1 ) 


= -G2(z)T^(z) 


(6-2) 
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•K^Gj^(s)Aj(s) .K^K^G^(s) 


Tc(z) 


-G| (2)T^(z) 


where 


(s) = 1 + K^s ‘ + K^s"^ 

A(s) == 1 - + K2S~’ + K^s"^ 


Gj(z) = 2,7|^ 

G,(z) . 2.78.<10-V(z^ I) 
2(z - ]f 


(6-3) 




(6-5) 


( 6 - 6 ) 


(6-7) 


A digital signal flow graph portraying Eqs. (6-l), (6-2), and (6-3) is 
shown in Figure 6-1. The characteristic equation of the system is written 


directly from Figure 6-1. 


K T 

A(z) = 1 + G 2 (z) + KjG^ (z) = 0 


(6-8) 


To obtain the equivalent transfer function that Q (z) sees, we divide both 


sides of Eq. (6-8) by the terms that do not contain Q. (z) . We have 

Go (z) K + ^ 

I . K,6(U) ' ° ■ 

Thus, ■ G.,(z) K + — 

r ( 7 \ = 2 I o z - 1 

eqo^ (z) 


C6-9) 


(6-10) 


Quantizer in the Channel 

Using the same method as described in the last section, let Q^ (z) denote 
the discrete describing function of the quantizer Q^ When only Q^ is 
considered effective, the following equations are written directly from Figure 


0^(z) = -G 2 (z)T^(z) 

= -Gj (z)T (z) 
r K. T • 


( 6 - 11 ) 

( 6 - 12 ) 


o z 


— Q^(z) + KjQ,(z) ^(z) 
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The digital signal flow graph for these equations is drawn as shown in 
Figure 6-2. The characteristic equation of the system is • 


A(z) = 1 + K^G] (z)C)j (z) + G^iz) 


K|T 


K + , 

o z - 1 


= 0 


Thus, the linear transfer function Q.(z) sees is 


(z) 


'"eql 


1 + G^iz) 


K + — 
o z 


^1^ 


(6-14) 


(6-15) 


CLuantizer in the T Channel 

— Q 


If Q..J. is the only quantizer in effect, the following equations are 
written from Figure 5~2. 

0^(z) = -G2 (z)Q.^.(z)T^(z) (6-16) 

51^(z) - -G^ (z)a^(z)T^(z) (6-17) 


K,T 


0^(z) + K,J2^(z) 


( 6 - 1 - 8 ) 


The digital signal flow graph for these equations is drawn in Figure 6-3. 
The characteristic equation of the system is 


A(z) = 1 + K^G^(z)Q.j.(z) + G 2 (z)ay(z) 


K|T , 


K + , 

o z - 1 


= 0 


(6-19) 


The linear transfer function seen by Qj(z) is 




K|T 


K + , 

o z - 1 


( 6 - 20 ) 


The discrete describing function of quantizers has been derived in a 
previous report^. By investigating the trajectories of the linear equivalent 
transfer function of Eqs. (6-l0), (6-15), and (6-20) against the critical 
regions of the discrete describing function of the quantizer, the possibility 
of self-sustained oscillations due to quantization in the digital IPS system 
may be determined. 



Sh 



Figure 6~3. Digital signal flow graph of IPS with Q.j in effect. 

Let denote the period of the self-sustained oscillation, and 

T = NT ■ 

c 

( 

where N is a positive integer ^2. T represents the sampling period in seconds. 

When N = 2, the periodic output of the quantizer can have an amplitude of 

kh, where k is a positive integer and h is the quantization level. The critical 

region of the quantizer for N = 2 is shown in Figure 6-^. For N = 3j the periodic 

output of the quantizer is a pulse train which can be described by the mode 

(•^o’ ’ ^- 2 }' and k^h are the magnitudes of the output pulses 

during one period. Similarly, for N = 4, the modes are described by (k^, k| ) . 

The critical regions of the quantizer for N = 3 and N = 4 are shown in Figures 

6-5 and 6-6, respectively. Figure 6-/ illustrates the frequency loci of G (z) , 

eqo 

*^eql^^^’ *^eqT^^^ (6-15), and (6-20), respectively, together 

with the corresponding critical region of the quantizer. The value of K| is 
equal to 10^ in this case. it so happens that the frequency loci of these 
transfer functions are almost identical. Notice that the frequency loci for 
the range of O.O 6 _< T £0.18 sec overlap with the critical region. Therefore, 
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sel f ~susta i ned oscillations of the mode N = 2 are possible for the sampling 
period range of 0.06 ^ T ^ 0.18 sec. It turns out that the frequency loci for 
N = 2 are not sensitive to the value of K| so that the same plots of Figure 
6-7 and the same conclusions apply to K| = 10^ and Kj= 10^. 

Figures 6 - 8 , 6 - 9 , and 6-10 illustrate the frequency loci of G (z) , G , (z) 

eqo eql 

6 ^qy(z) for N = 3 and for Kj = 10 , 10 , and 10 , respectively. From 
Figure 6-8 we notice that for Kj = 10^ the frequency loci do not intersect with 
the -critical region for any sampling period, .Thus, for K| = 10^ the N = 3 mode 
of oscillations cannot occur. 

For Kj = 10^, Figure 6-9 shows that sustained oscillations for N = 3 
would not occur for T < 0.075 sec and large values of T. For .Kj = 10^, 

Figure 6-10 shows that the critical val’ue of T is increased to approximately 

0.085 sec. 

For N = ^, Figures 6-11, 6-12, and 6 -I 3 illustrate the critical region 

c 6 7 

and the frequency loci for Kj = 10^, 10 , and 10 , respectively. In this case, 
self-sustained oscillations are absent for Kj = 10^ for any sampling period. 

For Kj = 10^, the critical sampling period is approximately 0.048' sec for 
quantizers Qj and whereas for 0^ the critical T is 0.0/ sec. The stability 
condition is improved when Kj is Increased to lO , for and 0..^., the critical 
values of T are 0.05 sec and 0.055 sec, respectively; for it is O.O 9 sec. 

As N increases, the critical regions shrinks toward the negative real 
axis of the complex plane, and at the same time the frequency loci move away 
from the negative real axis. Thus, the N = 2, 3, and 4 cases represent the 
worst possible conditions of self-sustained oscillations in the system. 

The conclusion of this analysis is that the sampling period of the IPS 
system should be less than 0.048 second, in order to avoid self-sustained 



oscillations excited by the quantizer nonlinearities described in these 
sections. 



Figure 6-^. Critical region of -1/Q.(z) of quantizer for N = 2. 
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7. Digital Computer Simulation of the Digital IPS Control System 

With Quantization 

The digital IPS control system with quantizers has been studied in 
Chapter 5 and 6 using the gross quantization error and the describing function 
methods. In this chapter the effect of quantization in the digital IPS is 
studied through digital computer simulation. The main purpose of the analysis 
is to support the results obtained in the last two chapters. 

The linear IPS control system with quantization is modeled by the block 
diagram of Figure 5~1, and it is not repeated here. The quantizers are assumed 
to be located in the T^, 0^, and channels. From the gross quantization error 
analysis it was concluded that the maximum error due to quantization at each 
of the locations is equal to the quantization level at the point, and it is 
not affected by the other two quantizers. In Chapter 6 it is found that the 
quantizer in the T^ channel seems to be the dominant one as far as self-sustained 
oscillations are concerned. It was also found that the self-sustained oscillations 
due to quantization may not occur for a sampling period or approximately 0.05 
sec or less. 

A large number of computer simulation runs were conducted with the quantizer 
located at the T^ channel. However, it was difficult to induce any periodic 
oscillation in the system due to quantization alone.. It appears that the signal 
at the output of the quantizer due to an arbitrary initial condition will 
eventually vary between \\^/l and -h.^/2 indefinitely in a random fashion. 

This points to the fact that the system the quantizer sees is not a low-pass 
filter so that the describing function method becomes inaccurate. However, 
the results still substantiates the results obtained in Chapter 5; i.e., the 
quantization error is +h^/2. Figure 1~\ and 7"2 show typical responses of 
the simulation runs for Kj = 10^ and T = 0.08 sec. When K = 10^ the system 
is unstable- 
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8. Modeling of the Continuous-Data IPS Control System With Wire Cable 
Torque and Flex Pivot Non 1 i near i t ies 

in this chapter the mathematical model of the IPS Control System is . 
investigated when the nonlinear characteristics of the torques caused by the 
wire cables and the friction at the flex pivot of the gimbal are considered. 

In Chapter 3 the IPS model includes the wire cable disturbance which is 
modeled as a nonlinear spring (Figure 3“3) . The combined effect of the wire 
cable and flex pivot is also modeled as a nonlinear spring characteristics as 
shown in Figure . 

i} 5 

In this chapter the Dahl model ’ is used to represent the ball bearing 
friction torque at the flex pivot of the gimbal, together with the nonlinear 
characteristics of the wire cables. 

Figure 8-1 shows the block diagram of the combined flex pivot and wire 
cables torque characteristics, where it is assumed that the disturbance torque 
at the flex pivot is described by the Dahl dry friction model. The combined 

'f 

torque is designated Tj^. 



Figure 8-1. Block diagram of combined nonlinear torque characteristics 
of flex pivot and wire cables of LST. 
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Figure 8-2 shows the simplified IPS control system with the nonlinear 
flex pivot and wire cable torque characteristics. 

The nonlinear spring torque characteristics of the wire cable are 
described by the following relations: 

^ Hy.^SGN(c) + (8-1) 


where is in N-m, in N-m/rad, e is in rad, and Ty^(e) in N-m. 

Equation (8-1) is also equivalent to 


T+ (e) - 

£ > 0 

( 8 - 2 ) 

“ ‘“WT 

e < 0 

( 8 - 3 ) 


It has been established that the solid rolling friction characteristic 
can be approximated by the nonlinear relation, 


where 


dTpp(e) 

de 


" ^FPO^ 


i = positive number 
Y = positive constant 


Tpp, = Tpp$GN(e) 


TppQ = saturation level of Tpp 


For i = 2, Equation (8-^) is integrated to give 


(8-4) 


e + Cj 


^fpo^ 

£ > 0 

( 8 - 5 ) 

r 

e + C 2 

T(Tpp + 

^FPO^ 

£ < 0 

( 8 - 6 ) 

are constants of 

! nteg rat ion , 

and 


T _ 

FP 

"^FP 


e ^ 0 

( 8 - 7 ) 

Tf“p = 

■'fp 


£ < 0 

( 8 - 8 ) 


The constants of integration are determined at the initial point where 



















n 


Ther 


= initial value of £ 

1 


T-d- = initial value of 
Fr I FP 


C, = -e. 


^^^FPi ~ "'^FPO^ 


£ > 0 


C. = 


-£. 

t 


"^^^FPi ^FPO^ 


E < 0 


(8-9) 

( 8 - 10 ) 


The main objective Is to investigate the behavior of the nonlinear elements 
under a sinusoidal excitation, so that the describing function analysis can be 
conducted. 


Let e(t) be described by a cosinusoidal function, 



e(t) = Acoswt 

(8-11) 

Then, 

e(t) = “Ao)sinwt 

(8-12) 

Thus, 

£. = -A £ > 0 

1 - — 

(8-13) 


£. = A ' £ < 0 

1 — 

(8-14) 

The constants of integration in Eqs. (8-9) and (8-1 O) become 



1 

' y(t - T T 

YUppj 

(8-15) 

C2 = ~ 

^ y(T ^~T ^ 

FPi FPO'^ 

( 8 - 16 ) 

Substitution of Eqs. 

(8-11) and (8-15) in Eq. (8-5) and simplifying, the 


solution of T|I^ is 

Vp R-1 

^FPO j(l - cosojt) + 

(8-17) 

• 

which is valid for e 

_> 0 or (2k+l) TT < tot < (2k+2)iT, k = 0, 1,2, . . 



a - 2yATppQ 

(8-18) 


P - 1 ,/a^ + 1 _ ■^FPJ 

(8-19) 

* 

Simi lar ly, for e < 0, 

using Eqs. (8-11) and (8-16) in Eq . (8-6), we have 
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^FP ^ 

"^FPO |(1 - coso)t) + ^ 


( 8 - 20 ) 


which is val i.d for 2kTr<^0Jt £(2k+l)Ti,k = 0, 1,2, ... 

' The express ions for and T^Tp obtained in Eqs. (8-17) and (8-20) together 
with those of 7 ^q(£) Eqs, (8-2)- and ( 8 - 3 ) are useful for the derivation of 

the describing function of the combined nonlinearity of the wire cable and 
the flex pivot characteristics. 

The torque disturbance due to the two nonlinearities is modeled by 


T+ = T + + T 
N WC FP 


= + Ky^AcosoJt + T 


+ |-( 1 - coswt) 

FPo a7I r~ 

j(l - coswt) + 


( 8 - 21 ) 


(2k+l)7r < (ot < (2k+2)TT, k = 0, 1 ,2, . . . 


T~ = T “ + T “ 
N WC FP 


- |-(1 - coswt) 

•‘^WT "^FPO IT ^ r 

-j{] - coswt) + ^ 


( 8 - 22 ) 


Figure 8-3 shows the Tpp/T^p^ versus e/A characteristics for several 
values of A when the input is the cosinusoidal function of Eq. (8-11). 
Figure 8-4 shows the normalized (3Tpp + / (STpp^ + + H^.p) versus 

e/A for several typical combinations of A, H^,p, and Ky.p. 








Figure 8-^. Normalized flex 
for IPS with cos 
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9. Describing Function of the Combined Wire Cable and Flex Pivot Nonl ineari t ies 


Figure 8-1 shows that the disturbance torques due to the wire cable and 
the gimbal flex pivot are additive. Thus, 


^WC ^ ^FP 


(9-1). 


For the cosinusoidal input of Eq. (8-11), let the describing function of 
the wire cable nonlinearity be designated as N (A) and that of the flex pivot 

W 

nonlinearity be Nj_p(A). Then, in the frequency domain, the total disturbance 
torque is 


T|^(o)) = Npp(A)e(w) + Ny^(A)e(co) 
+ (Npp(A) + N„^(A))£(w) 


WC 


(9-2) 


Thus, let N(A) be the describing function of the combined flex pivot and wire 
cable nonlinear characteristics, 

N(A) = Npp(A) + N^^(A) (9-3) 

The describing function of the Dahl sol_Ld_friction nonlinearity has been 
derived elsewhere for the cosinusoidal input. The results is 


Npp(A) = 


^1 - 

A 


where 


A, = - ^ T 


IT FPO ttAy 
C 


In 


C^ + A^ 


- A 


^1 


—=L 

/c j - 


- 1 


(9-M 

(9-5). 

(9-6)- 


The describing function of the wire cable nonlinearity Is derived as follows. 

For a cosinusoidal input the input-output waveform relations are shown in 
Figure 9-1* The wire cable torque due to the cosinusoidal input, over one 
period, is 

\c^^^ ^ KyjAcoswt - 


0 < wt < TT 
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Figure 9~1 . 


Input-output characteristics of wire cable torque 
non] ineari ty . 
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" Ky^Acoswt + 


7T < Wt < 2TT 


The fundamental component of the Fourier series representation of T,,^(t) 

WC 

Twc(t) = A^slnwt + B^cosojt 

=a/ A^ + cos(wt - <{)) (9-8' 

^ = tan — (9-9, 

^1 


1 

~ ~ T. ,„(t)slnwt dwt 

Jn 




The coefficients A^ and are derived as follows; 


2 

= — T, ,-(t)sInO)t du)t 
-IT Jq WC 

2 

~ ¥ j ~ Hy^)sInojt d(ot 


“ T (t)coso)t dcot 

IT Jq WC 


■ ^ C 


„Acos(0t - H, ,_)coswt dwt 
1 W I 


Then, 


1 1 
^ 1+^1 




-1 r^^wT 
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The describing function of the wire cable nonlinearity is written 


E - JA /aJ + bJ 


’WT' 


A 


A 








WT 




2 / -1 — yj 

+ /tan TTKy^A 


(9- 16') 


For the combined nonlinearity, the describing function is the sum of the 
two describing functions. However, since there are three ball bearings on the 
flex pivot,- the final expression is ' 

N(A) = N^^y(A) + 3Npp(A) = N^(A) + jN,(A) 

where 


»,(A) = 3 ^ 
ya 


N I (A) =3 


ttA FPO 


ttA y 


In 



(A) + 

1 


A^ 

“ 1 

r' 

+ A-|' 


- aJ 


I 


(9-17) 


^T 


4h 


WT 


TTA 


(9-18) 


(9-19) 


Asymptotic Behavior of -1/N(A) for Very Small Values of A 
The asymptotic behavior of -1/N(A) for very small values of A can be 
derived analytically. It can be shown that' 

.. linN^(A) (9-20) 


and 


Therefore, 


1 im N j (a) = 1 im 
A -^0 


4H, 


WT 


A-^0 


1 im 


TTA 


(9-21) 


-1 


A-.0 “ A-lo 


-1 


iH 


j 1 im 


WT 


= 0/-270 


(9-22) 


ttA 
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Asymptotic Behavior of -]/N(A) for Very Large Values of A 
For very large values ,of A it can be shown that 
lim N (A) - 

1 im N| (A) - 0 

A-tCO 


Then , 


1 im -1/N(A) = 


(9-23) 

(9-25) 


Magnitude Versus Phase Plots of -1 /N(a) of the Combined Non 1 inear j t ies 
A digital computer program for the computation of N(a) and -1 /N(a) is 
listed in Table 9“1 - The constant A is designated as E in this program. 

The parameters of the honl inearities are: 

TppQ = 0.00225 N-m 

Y = 3.2kkh ^ 10^ (N-m-rad)“' 

= 0.25 to 100 N-m/rad 
H, = 0.01 to 1 N-m 

Figure 5-2 shows the plots of -1/N(A) for the non 11 near! ties in magnitude 
(db) versus phase for varies combinations of K^.^. and It is seen that 

varying the value of between the range of 0.01 to 1 does not affect the 
curves appreciably. 


Prediction of Self-Sustained Oscillations in the IPS System With the ' 
Combined Nonlinearity By Means of the Describing Function Method 

The characteristic equation of the nonlinear IPS control system with the 
r 

wire cable Dahl-model nonl inearities is determined from Figure 8.-2. 

1 + N(A)G^^(s) = 0 
eq 

-where N (A) is defined in Eq. (9“17)j and 


(9-26) 



83 


TABLE 9-1 

IPS CONTINUOUS DESCRIBING FUNCTION OF COMBINED NONLINEARITY 


CflLCULHTIOH FOP -1-N 

comple:-: gvv jONjov 

FlEfiL+y P'.cIU ' jPI j PRD j TO j bflMMfl 'ESThPT jE 'R fl jP 'TGFI .’TGFN-' 
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G (s) - Q-0013946(s^ + Q.0012528s^ + 0. 0036846s) , 

eq ^ A (7) ~~ ’(9-27) 

where A-(s) is defined in Eq. (3-k) . I^tice that G (s) in Eq, (9-2.7) is equal 

a eq 

to S' times the G (s) in Eq, (3~3) , since in Chapter 3 the^ input to N is 

eq 

e whereas now It is e. The frequency, plots of G (s) of Eq. (9-27) are plotted 

in Figure 9*-2 for Kj = 10 and 10 . Similar to the curves In Figure 3“6, these 

frequency loci for G^^(s) have two equilibrium points for each curve, one stable 

and the other unstable. For instance, for K, = 10^, the G (s) curve intersects 

1 eq 

the -T/N{a) loci, at w = 0.138 rad/sec and oj = 0.055 rad/sec. The equilibrium 
point that corresponds to u = 0.138 rad/sec is a stable equilibrium point, 
whereas o) = 0.055 rad/sec represents an unstable equilibrium point. These 
results are very close to those obta i ned in Chapter 3 where the flex pivot 
is presented as a spring. Therefore, the impact of using the Dahl solid friction 


model is not great although all the loci are substantially different. For' 

K|= 10 , the stable equilibrium point is at o) = 0.16 rad/sec, and the unstable 
equilibrium point is at tu = 0.05 rad/sec. 

Figure 9-2 shows that for the system paramters used, the intersections 
between G^^(s) and -1/N(A) all fall on the portion of the -1/N(A) loci that lie 
on the -270° axis. This means that as we vary the values of and of 
the wire cable nonlinearity characteristics within the stipulated ranges, only' 
the amplitude of oscillation. A, will be varied. Equations (9-21) and (9-22) 
further show that for small values o^ -1/N(A), which correspond to the range 
of intersections in the present case, the amplitude of oscillation is not 
sensitive to the values of y, T^pQ, and K^.^. However, the amplitude, of 
oscillation. A, is directly proportional to the value of Typical results 

of the stable sustained oscillations are tabulated in Table 9"2. 
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Figure 9-2. Frequency response plots and describing function, loci of IPS system with 
flex pivot and wire cable non 1 i nea r i 1 1 es (Dahl model). 










TABLE 9-2 


s 

V 

u 

WT 

A (rad) 

(arc-sec) 

01 (rad/sec) 

10^ 

0.25 

0.01 

10-^ 

0.2 

0.138 

105 

0.25 

0. 1 

10-5 

2,0 

0.138 

10^ 

1.00 

0.01 

lo"^ 

0.2 

0.138 

10^ 

1.00 

1 .'00 

10“^ 

20.0 

0.138 

lo5 ' 

5.00 

0.01 

lo“^ 

0.2 

0.138 

105 

5.00 

1.00 

lo"^- 

20.0 

0.138 

105 

10.00 

0.01 

lo"^ 

0.2 

.0.138 

105 

10.00 

1.00 

io“^ 

20.0 

0.138 

lo5 

25.00 

1.00 

10"'^ 

20.0 

0.138 

10^ . 

.1 00.00 

l.OO 

lo"^ 

20.0 

0.138 

10^ 

0.25 

0.01 

3xlo"^ 

0.006 

0.16 

10^ 

0.25 

0.1 

3xl'0"^ 

0.06 

0.16 

. 10^ 

1 .00 

0.01 

3x1 o"^ 

0.006 

0.16 

I0<^ 

i.oo 

1.00 

3xlO"^ 

0.6 

0.16 

10^ 

5.00 

0.01 

3xl0"® 

0.006 

0.16 ■ 

10^ 

5.00 

1.00 

3x10“^ 

0.6 

0.16 

10^ 

10.00 

0.01 

3xlo"^ 

0.006 

0.16 

io^ 

10.00 

1.00 

3x1.0“^ 

0.6 

0.16 

• 10^ ■ 

25.00 

1.00 

3x10"^ 

0. 6 

0.16 

10^ 

100.00 

1.00 

3x10“^ 

0 . 6 

0.16 
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with those obtained 'in 

30, the following comparisons 



Results in 

Chapter 3 

Dahl model 

resul ts 


A (rad) 

w (rad/sec) 

A (rad) 

(jd (rad/sec) 

lo" 

I 

3 . 17 x 10 '^ 

-0. 1 6 

3x10"^ 

0.16 

10^ 

8.19x10~^ 

0.14 

0 

1 

0.138 


It is of interest to compare these results 

I 

Chapter 3- Using the results tabulated on page 

I 

are obtained. 

I 

= 1, Ky.|. = arbitrary 


Therefore, we see that for all practical purposes these results are 


i den ti cal 
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10. Modeling of the Solid Rolling Friction by the First-Order Dahl-Model 
it has been established that the solid rolling friction characteris- 
tics can be approximated by the nonlinear relation 

de “ ^^^FPI ” ^FPO^ OO-l) 

where i = positive number 
y = constant 

Tpp I - TppSGN(e) 

Tpp(e) = friction torque 

TppQ = saturation level of Tpp 
e = angular displacement 

The describing function of the friction nonlinearity for n = 2 has 
been derived. In this chapter the input-output relationship will be ob- 
tained by solving Eq. (10-1) with i = 1. The describing function for the ' 
i = 1 case is deriv^ed in the next chapter. 

Let the ‘angular displacement e be a cosinusoidal function, 

eit) = A cos wt . (10-2) 

Then, we can write 

dTpp(e) dTpp.(e) . 

dt ^ d£ ^Tpp, - Tppg) (10-3) 

where i has been set to 1 . 

Since Tppi =TppSGN(e), Eq.(l0-3) is wr i tten 

dTpp(e) 

^ -yAo) sin tot (Tpp - Tppg) e > 0 

= yAw sin tot (Tpp -f- Tpp^) e < 0 (10-4) 

For e-< 0, 2 ttI< < wt £ (2k -t- 1 )tt, k = 0, 1 ,2, . . , , 
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Equation (10-4) is integrated on both sides to give 


T [c(t)] , f . 

tP dTpp(e) /CO 

Tpp[e(0)] ^FPO^ Jo 


fCOt 

= yAu) s i n (OTdx 
■ n 


Thus , 


T [E(t)] (jOt 

^ri(T_p + T ) = -yA cos cot 

Tpp[dO)] 


£n(Tpp(e) + TppJ - )ln(Tpp(0) + T^pJ = -yA(cos cot - l) (10-7) 


Tpp(£) + TppQ 

Tpp(O) -I- Tppj^ 


= -yA(cos cot - 1) 


(10-8) 


Since for the cosinusoidal input e(0) = A and e _< 0 for 27rk ^ cot < 
(2k + l)'n‘, k = 0,1,2,..., T^p (O) = T^p | ^ 0 . This is because at £ = 0 , 
£ is decreasing, and Tpp acts in the direction opposite to the motion; 
thus ^FPI - Equation (10-8) is written as 


Tpp ( e) + Tpp„ 

An — T= = -yA(cos cot - 1) 

FPI FPO 


(10-9) 


. PPI I , -yA(coscot-l) 
T T ^ > e 

FPO 'fPO 


( 10 - 10 ) 


For £ > 0, (2k + 1 )tt < cot < (2k + 2)-fr, k = 0, 1 ,2, . .. . 
Equation (10-4) is integrated on both sides to give 


Tpp(2ir) 


dT^p(e) 


Tpp[£(t)] FP " 'fPO^ Jcot 


-yA s i n COT dcoT 


(10-n) 
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Carrying out the integration, we get 

= “YA(1 - cos (ot) 


£n 


Tpp^e) ■ Tj.pQ 




( 10 - 12 ) 


Now at tot = 2 tt, T|_p(27r) = Tj_p ^ ^ 0. Equation (10-12) leads to 


T,p(s) 


T 


FPI 


— 1 + 


FPO 


T 


FPI 


[^FPO 


yA(cos (jJt - 1 ) 


(10-13) 


Let 


"‘Vp/'^fpo 

Then, Eqs. (10-10) and (10-13) become 


T fr) 

= -1 + (R + 

FPO 


(10-11{) 




FP YA(costot-l) 


= 1 + (r - 1), 


(10-15) 


( 10 - 16 ) 


FPO 


respectively. 

In order to evaluate R, we equate the last two equations at cot = IT. 

Then, 

-1 + (R + De^"^^ = ] + (R - l)e"^^^ 


(10-17) 


The solution is 


R = 


„ 2yA -2yA 

2 e ' + e 

^2yA _ g"^YA ^2yA _ g“2yA 


or 


R = csch (2yA) - cot h (2yA) 


(lo-rs) 


Since e = A cos wt and e. = e(0) = A, Eqs. (10-15) and (10-16) are written 


as 


Tpp(e) 

''^FPO 


+ (R + 1 ) e 


-Y(£-£. ) 


e < 0 


(10-19) 
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y(£“£.) 

= 1 + (R - Oe ' 


e > 0 


Figure 10 1 shows the T^p/T^pQ versus e/A characteristics for several 


values. of A when the input is the cosinusoidal function of Eq. (10-2). 
Figure 10-2 shows the normalized (3Tpp + versus 

e/A for several typical combinations of A, and 
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11. Describing Function of the Fi.rst-Order Dahl Model Solid Rolling' 
Friction 

The mathematical description of the first-order Dahl model of the 
solid rolling friction is presented in the last chapter. The frictional, 
torques for the two ranges of £ for a cos i.nusoidal input displacement are 
given by Eqs. (10-19) and (l0-20). These equations are rewritten in the 
fol lowing form: 

Tpp(.) = - I) (in-,) 

TppOd = - ,) ' (,,-2) 

For the cosinusoidal input, s(t) = A cos (ot , let T^_(e) be approxi- 

rr 

mated by the fundamental component of its Fourier series, representation, 
i . e- , 

Tpp(e) = sin ojt + cos tot (11-3) 

The describing function of the friction nonlinearity (i = 1) is defined as 

B] - jA 

Nfp(a)=-^-^^ ^ (11-4) 

where 

1 

A, = ■— T_„(e) sin tot dwt 

I TT J l-F 


1 

TT 




((Tppi + T^pq) 


-'YA(costot-l ) 


- T^pq) sin o)t d(jJt 


rZTT 


+ — 
7T 


' 7T 




YA(co 50 Jt-l) 


+ T^pq) sin tot doot (11-5) 


Evaluating the integrals in the last equation, we have 


-4t 




±P0 + _]_ fy (e"2YA 2yA _ ^) + T (e^^^ - 
IT ttyA [ FPr ' FPO^ ' 


(U-6) 
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or 


“4T r 

A = FPO ^ _2_ 

1 TT fryA 


Tpp|(cosh 2yA - 1) + Tpp^(sinh’ 2yA) 


(11-7) 


1 


" f(T„, + - T 


®1 TT Iq ';^''FPI ■ ’FPO 


FPO 


4 - — 
IT 


2ir 


TT 


(T - T )gYA(coso)t-l) 

^ FPl FPO^ 'fPO 


In order to evaluate the i ntegra 1 s. of let us represent e 
power series: 


cos o)t dwt 

cos o)t dojt ( 1 1 -8) 
yAcoswt 


as a 


B1 


-YAcoswt , 

e = 1 - yA cos wt 


Consider the integral 
■IT 
■ 0 


2 2 ^ ^ 

(yA) cos (ot _ (yA) ^ cos^ cot 

2 1 ***^ ^ I • 


3! 


-yAcosojt 

e cos 03t dwt 


"" ' (yA)^ 2 

1 - yA cos (ot + , cos wt 


0 


2! 


Mil eos3 + 


3! 


S i nee 


•ir 


cos tot dwt = 0 for m = odd integers 


(n-9) 


] cos (Ot doJt 
! ( 11 - 10 ) 


( 11 - 11 ) 


Eq. (11- 10) becomes 


I — j cos*^* (Ot d(0t i = odd integers (11-12) 


I rTT 


B1 


i=0 


I ! 


i + 1 


0 


Evaluating the integral, the result is 
, _ -iryA 


B1 


(yA)^ 

1 

(YA)'* 


3 

, (yA)^ 

f 1 -3 -51 

2! 

l^J 

if! 

IsJ 


6! 

_ 8 • 6 ♦ if 


+ 

8! 


8 


1 ♦ 3 • 5 • 7 
10 • 8 • 6 • A 


+ . . 


(11-13) 
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SimHarly, the integral 


B2 


2tt 


TT 


yAcoswt ^ , 

e cos tot dtot 


is evaluated, and the result is 


*B2 "'bI 


Ol-U) 


Substituting the results of 1^^ and 1^^ into Eq. (11-8), we have 


B, = 


^FPI ^FPO yA , ^FPI " ^FPO -yA, 

TT 'b1 TT ^ 'B2 


- Ill [t «"TAi , / YA , -yA\ 


( 11 - 15 ) 


For very small values of yA 


1 ~ ^ya 

Bl 2 


( 11 - 16 ) 


Equation (11-15) becomes 




( 11 - 17 ) 


or 


B, £ -yA 


Tppi sinh (yA) + TppQ cosh (yA) 


(11-18) 


For large values of yA, 1^^ becomes very large. However, we shall 
show that Tpp I approaches “fppQ as yA becomes very large, so that B^ 
becomes zero. 

We shall now investigate the limiting values of A^/A and B^/A when 
A approaches zero and infinity. Since 


1 im T p, = 0- , 
A^O 


( 11 - 19 ) 


I I m 
A -^0 


1 i m 
A-^0 


"^^FPO"^ ^^FPO"^ 

“T 


2ttAy 2TrAy 


= 0 


( 11 - 20 ) 
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1 im ~ = 
A^O '' 


( 11 - 21 ) 


Therefo re, 


1 im -1/N (A) 

A->0 ^ ' 


When A approaches infinity. 


yt 


FPO 


(f]-22) 


l!Z ' '"fpo 


B, 


lim ^ = - 1 

A-+~ 

= 0 


-T e^^ + T 
FPO FPO 


B1 


(11-23) 


(n-2if) 


The value of A^/A also approaches zero as A becomes very large; 
however it decreases at a much slower rate than B^/A. Thus, 


1 im 

A-wo 


-l/N,p(A)] = ^ 


Figure 11-1 shows the plot of -l/N|_p(A) in the magnitude (db) versus 
phase coordinates for 


Y = 13429,75 (N-m-rad)"^ 


TppQ = 0,0088 N-m 


Magnitude versus Phase Plots of -l/Npp(A) of the Combined Nonlinearities 

■For the combined nonlinearity of the Dahl friction model (i = 1) and 
the wire cable, the nonlinear describing function is written as (Eq. (10-17)) 


N(A) = N|^^(A) + 3Npp(A) 


(11-25) 


where N^j(A) is given in Eq. (9-16). 

Figure'|fl-2 shows the plots of -1/N(A) for the combined nonl.inearity 

in magnitude versus phase for various combinations of K,,_.and H 
. / wT 


wr 


These 



CIBELS -1/N(A) 
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'n(A) 


K,=10 


I'V,. " 100, 

a) V “ 50- = I 

(3) = 25, “ 1 

(“I) l\jj - 10, = 0.01 

(5) Ky^ = 5, Hyj. = 0.01 

^/T “ 5. Hy^ » 0. 1 

(7) ><KT “ >^„T “ 

(8) Ky^ = 0,25,Hy^ = O.I 



-oj roaTse^ 

Prio’io 


lo'lo^ -1 




/0.07 . 


32^' -300“ t 280" -260° -240“ -22?^20d’ _14Cf -l^d* -lOtf "" -gtf -6tf '-4C? 

Figure, n-2. Frequency response plots and describing function loci of IPS system with 
flex pivot and wire cable no'nl 1 near! t ies (Dahl model i = 1). 
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curves are similar to the plots shown in Fig. 9-2 which are for i = 2 in 
the Dahl model, especially when the values of A are very small and very 
large. 

The frequency loci of G (s) of Eq. (9~27) are plotted in Fig. 11-2 

0CJ 

5 6 

for Kj = 10 and 10 . Similar to the cases in Fig. 9“2, these frequency 
loci have two equilibrium points for each curve, one stable and the other 
unstable. 

For K| = 10 the frequency of oscillation at the stable equilibrium • 
is approximately 0.16 rad/sec, and is rather independent on the values of 
Kyy and This result Is identical to that obtained in Chapter 9 when 

i = 2 is used for the Dahl friction model. 

C 

Figure 11-2 shows that for Kj ~ 10 the frequency of oscillation at 

the stable equilibrium varies as a function of the values of K^^y and H^y. 

For the various combinations of K^y and H^y shown in Figure 11-2, the 

variation of frequency is not large, from 0.138 to 0.1^ rad/sec. Of more 

importance is perhaps the fact that when i = 1 , the amplitude of oscillation 

A is larger as compared with that for i = 2. For example, for K^y = 100, 

^WT ~ Kj = 10 , Figure 9~2 shows that the amplitude of A is approximately 

“5 

10 for k = 2, whereas for the same set of parameters Figure 11-2 shows 
that A = lO"^ for i = 1 . 
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12. Modeling of the Solid Rolling Friction by the ith-order Dahh Model 
and The Describing Function 

In the previous sections the solid rolling friction was modelled by Eq. 
(10-1) with i = 1 and 2, In general, the exponent i can be of any other value. 
In this section we shall derive the mathematical model of the solid rolling 
friction for i 1 . - . 

Let the frictional characteristics be approximated by the nonlinear 
relation 


dTpp(£) 

de “ 


( 12 - 1 ) 


where 


i = positive number ^ 1 
Y = constant 
Tpp, - Yp SGN(s) 

Tpp(e) = friction torque 

TppQ = saturation level of Tpp 
e = angular displacement 

Let the angular displacement e be a cosinusoidal function, 
e(t) = A cos 0 )t 

Then, e(t) = -Aw sin wt 

We can write 


( 12 - 2 ) 


dT 


FP 


dt * - - W 

In view of Eq. (12-2), the last equation is written 

dT. 


(1-2-3) 

(12-4) 

(12-5) 


FP 


dt 


= - yAw sin o)t(Tpp - T^p^) 

- - yAw sin ^t(-T^p - T^p^) 


e ^ 0 
£ < 0 


If e ^ 0 for 2k7r - 5 ut < (2k+l)TT , k = 0, 1,2, '. . 


(12-6) 

Eq. (12-6) is integrated 


-'dT 


FP 


Tpp(e(0)) 


(t - t ) ' 

^ FP FPO^ 


= - yA 


■wt 


sin udu = - yA(~cos u) 
ojt = 0 


wt 


REPRODUClBir^|TT^«|' 'TltS) 
ORIGffNiS. PAGE « POOBi' 
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= - yA(~cos ut + ]) 


( 12 - 7 ) 


The last equation becomes 


(T - T ) 
^ FP 'fpO^ 


-i + 1 


[ (- i+1 ) 


FP (T - T ) 
^ FP FPO^ 


-(i-1) 


^pp j ”(■“!) 


(T - T ) 
^ FP i FPO--^ 

mi 


(i-1) 


= YA(cos wt “ 1 ) 

where T^p. = Tpp(e(0)) and T^p = Tpp(e(t)). 

Equation (12-8) is further simplJfi.ed to 

^^FP ” ^FPO^ ' ~ ^^FPI ~ ^FPO^ ^ ~ -(i-l)YA(cos o)t - 1) 

Defining 

Eq. (12-9) leads to 

T. 


(12-8)- 


= ^FPi^^PO 


■Cl 2-9) 
(12-10) 


FP 

FPO 


= 1 + 


R - 1 


{1 “ (i-l)YA(TppQR-l)}^‘‘^^cos wt - 


( 12 - 11 ) 


if £ < g, for (2k+l)Tr < 0 )t < (2I<+2 )tt, k = 0, 1, 2, . . . ; Eq. (12-6) Ls 

integrated to give 
(Tpp [e(27r)) 


dT 


FP 


Tpp(s(tot)) 


(-T - T ) ' 

^ FP FPO^ 


2ir 


YA sinu du = A(T - cosut) 


( 12 - 12 ) 


cot 


Following the same steps as in Eqs. (12-8) through (12-11), we have . 

U + 1) ... ,Ml/(i-l) ^ (12-13) 


FP 


R + 1 


FPO {1 + (i-l)YA(-TppQ(R+l))’'’‘"'^ (cos wt - 1 ) } 

In order to evaluate R, we equate Eqs. (12-11) and (12-13) at cot =ir. After 
simplification, the result is 


2 + R 


(1 + S(R-l) (l-l)jl/('-l) {!.- a'[-(R+l)J 


[0 + S(R-1)^''?V^^‘“^^ {1 - S[-(R+1))^'"^^}’/^''’^ 

(i-1) 


(I2~li}) 


a = 2(i - DyAT 


FPO 


where 


(12-15) 
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Once the value of i {\ ^ 1) is' sped f led , R can be solved from 
We can show that when i = 2, 

S - a - 

which is identical to Eq . (8-18), and 




which is the same result as in Eq. (8-19). 

Once R is determined from Eq. (12-14), the torque relationships 
expressed by Eqs. (12-11) for £ > 0 and Eq. (12-13) for e < 0. 


Describing Function For the ^Dahl Model For i ^ 1 
Let 


$ = */2 = (! - DtAT^^o'' ' 
Equations (12-11) and (12-13) are simplified to 


R 


FP 

(1 - B(R - 1)^'“^^ (cos mt - 


FP 


R + 1 


+ 1 (e > 0) 


- 1 (e < 0] 


^■FPO {7 + 3[-(R + l)]^'‘’^{co5 tot - 1)}-’/^'-’^ 

For the cosinusoidal, input of Eq , (12-3), let the output torque 
represented by the fundamental components of its Fourier series; i.e, 


where 


and 


^pp " sin tot + Bj cos cot 
r27T 


N =- 

1 7T 


1 


0 

2tt 


Tpp sin cot dtot 


= - . T|_p cos cot dcof 

Substitution of Eqs. (12-19) and (12-20) into Eq; (12-22), we get 


A, . I™, f f R + I 

1 IT J 


oi{l + e(-(R + l))^‘"’^cos cot - ' 


Eq. (12-14). 
(12-16) 

(12-17) 

are 

(,12-18> 

(12-19) 

« 

( 12 - 20 ) 

be 

( 12 - 21 ) 

(12-22) 

(12-23) 


IT 


si n cot dtot 



10 ^ 


Tppo r^TTr 


IT 


TT 


R - 1 


[(] - 6 (R - 1 ) 0 )t - D) 


+ 1 


si n (jot dcot 


The iast equation is -reduced to the following form: 


„ ^^FPO . 

A I — + 

1 7 T 7 T 


s i n wt dut 


0 {l +B(-(R + l)] (cos wt - l)..}'^ * ^ 


TppoCR - 1.) r2ir 


ir 


Let 


and 


Then 


sin wt dwt 

■FT [l + (R - 1 ) ^ ' ^ (cos wt •- 1 )]^^^' ^ 


x = 3(R'~ 0^' cos wt 


y= 3 (“(r + 1 ))^' '^cos wt 


dx = -B(R - 1 ) ^ ^sin- wt dwt 

dy + B(R + 1 }^' ^^sin wt dwt 
Equation ( 12 - 24 ) is written 

IjJ T j. i\ f3(R+l)^' 


^ = - 


FPO ^ Tfpq(R + 1 ) 

“n" _r, r n ! 1 \ ( i ~ 1 ) 


tt3(R + 1) 


. dy 


-3(R+1)^‘"’^ {1 +3(“(R+Oy^' + y} 


Tfpq(R - 0 


r 3 (R-l) 


it3(R - 1) 


TmT 


dx 


- 3 (R- 1 ) 


(i- 1 ) (l - 3 '(R- 1 )^' 


Let 


9 = 1 + 3 (- (r + 1 ^ ^ + y 


5 ^ = 1 + 3 (R - 1 ) 


(i- 1 ) 


Then m. d 9 = dy 

d 5 ^ = dx 

Substitution of the last four equations into Eq .• ( 12 - 29 ) yields 


A = - ^FPO^*^ 

TT 


t-3[~(R+0] J 


l+23(R+l) 


d9 


(9) 


l/(i-l) 


( 12 - 23 ) 


( 12 - 24 ) 

( 12 - 25 ) 

( 12 - 26 ) 

( 12 - 27 ) 

( 12 - 28 ) 

I 

TTrrry 

(12-29) 

(12-30) 

(- 12 - 31 ) 

( 12 - 32 ) 

( 12 - 33 ) 
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7Te(R “ 1) 


TmTJ 


1-26(r-1) 


(i-1) 




(^) 


l/(i-l) 


02-3^i) ‘ 


The ihtegralson the right-hand side of the last equation are now- carried 
out, and after s impl i float ion , -the result is 


A, = - 


FPO 


IT 


Tppo(R + 1) 

f(l + 26f-(R+l)l 

j(i-l)|(i-2)/(i-l) _ 

it3(r + 1) 
"’'fpo^^ ~ 

(i - 2] 
- 2f(R - Dl' 

i/U - 1) ■ J 

ttB(R - 1) 

(I - 2)/(i - 1) ■ 


(12-35) 


The Fourier coefficient B, Is determined as follows: 


■ 


’’fpo 


TT 


FPO 

ir 


0 

IT 


R + 1 


{l + (-B(R + 1)) ' ^ (cos wt - ])} 
R - 1 


i/(i-i) ; ' 


cos fait do)t 


{1 - e(B - 1) (cos mt - 


+ 1 


cos wt dwt ( 12 - 36 ) 


The last equation is simplified to 


T 


^1~- 


FPO 

IT 


f ' (R + 1) 


fir 


gf"(R + ])1 ^'~^^cos(ot dcot 


6[-(R+l)}^'~^^ Jo jT- 3[-(R + 1)}^'"’^ + 3(-(R + V)]^‘“’^coso)t}’^^' 


Tj 


(R - 1) 


3'(r'- 1) 


(T^ 


fl-n 


'IT 


B(R ~ 


1 ) ^ ‘ '^cos 03t dwt 


(1 + B(R - 1) 




- 3(R 


D^'-’^cos 


wt] 


1/Ti-l) 


(12-37) 


Letting 

X = B[-(R + 1)}^' ^^cos wt 
y = B(R “ 1 )^' ^^cos wt 
dx = -6(-(R + 1)]^' '^sin wt dwt 
dy = -B(R - 1 ) ^' '^sin wt dwt 


( 12 - 38 ) - 
(12-39) 
(l2-i<0) 
( 12 -^ 1 ) 


Eq, ( 12 - 37 ) becomes 
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B = f -(R + 1) 

1 •rr 


x(7t) 


cot ut dx 


&(-(R + D) ^x{0) {l - 6(-(R + 0) 


R 


ry(2Tr) 


3(R - 1) 


TT^ 


cot tilt dy 

U-1) 


yW [l + 3(R - I) y) 


For the first integral in the last equation, 


_ (12-42) 


cot (jOt 


and for the second integral, 

y 


cot wt = 


- I) 




Therefore, 
T, 


B, = 


FPO 

IT 


-(R + 1) 


fX ( it ) 


e[-(R + D] 


rrrr 


xdx 


(12-43) 


(12-44) 


(R - 1) 


e(R - 1) 


TT^ 


y(27T) 


Jx(0) {3^(-(R+l))^^'“’^-x^}’^^{l-3(-(R+l)}^”'^^+x}’^^' 




y(ir) (|3 (R-l) ) '■^(l + g(R-l ) ' ' ‘ ' -y) 


These integrals can be carried out only if the value of i is given (i ^ l). 
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13- Digital Computer Simulation of the Continuous-Data Nonlinear IPS 

Control System With Dahl Model 

The IPS control system with the nonlinear flex pivot torque modelled by 
the Dahl solid friction model is simulated on the digital computer for i = 1 
and i = 2. The block diagram of the IPS system is shown in Fig. 8-2, and' the 
nonlinearities are modelled by Fig. 8-1. 

The main objective of the computer simulation is to verify the results on 
the sustained- osci 1 lation predicted- by the describing function method. > 

Dahl Model i = 1 

The computer program using the IBM 380 CSMP for i = 1 in the Dahl model 
is given in Table 13~1« The simulation runs were able to predict and verify 
the results obtained by the describing function methed of Chapter 9. The 
difficulty with the long response time of the iPS system still exists in this 
case. Generally, it would be very time consuming and expensive to wait for , 
the transient to settle completely in a digital computer simulation. Figure 
13-1 shows the response of e(t) over a one-hundred second time interval, with 
e(0) = lO"^, e(0) = 0, Kj = 10^, = 100, = 1. For the Dahl model, i = 1, 

TppQ = 0.0088 N-m, and y = 13^29.75. The parameters of the linear portion of 
the system are tabulated on page If in Chapterl . Figure 13~1 shows that the 
response is oscillatory with an Increasing amplitude, and the period is 46 sec • 
or 0.136 rad/sec. Since it would take a long time for the amplitude to settle 
to a final steady-state value, we selected another initial value e(0) and 
repeated the simulation. Figure 13~2 shows the response of £(t) with e(0) = 

-3 

10 which is decreasing in amplitude as time Increases. Therefore, the stable 

-3 -5 

operating point should be at an amplitude between 10 and 10 , and the freq- 

uency of oscillation is 0.136 rad/sec. In general it would be very difficult 
to find the initial state which corresponds to the steady-state oscillation 
exactly. The results predicted by Fig. 9-2 are very close for the amplitude 
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Figure 13“2. Tine response of e(t) for 'Kj = 10^, = 100, 

e(0) = lO"^, i = 1, TppQ = 0.0088 N-m, y = 13^29.75. 
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Ill 

and 0) = 0.138 rad/sec. 

Dahl Model i = 2 

Table 13“2 gives the computer simulation program for the i = 2 case, 

with TppQ = 0.00225 N-m and y =" S2hhk.. All other system parameters are the 

same -as the i = 1 case. Figure 13“3 shows a stable response when the initial 

**10 

state e(0) is small, 10 . As shown in Fig. 11-2, when the initial state is 

small the stable equilibrium point is e = 0. Figure 13-4 shows another stable 

_o 

response which would take longer time’ to die out, when e(0) = i0~ . Figures 
13“5 and 13-8 show a sustained oscillation solution with apl'itude lying between 
10 ^ and 7 ^10 and a period of 44 sec or 14.28 rad/sec. These results are 
again very close to those predicted in Fig. 11-2. The simulations for the 
i = 2 case are carried out with K^.^. = 0.25, = 0.1, and Kj = 10^. 
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Figure 13~A. Stable response of E(t) for 
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